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Abstract. We present the first diffraction-limited speckle 
masking observations of the oxygen-rich AGB star AFGL 
2290. The speckle interferograms were recorded with the Rus- 
sian 6 m SAO telescope. At the wavelength 2.11 /j,m a res- 
olution of 75 milli-arcsec (mas) was obtained. The recon- 
structed diffraction-limited image reveals that the circumstel- 
lar dust shell (CDS) of AFGL 2290 is at least slightly non- 
spherical. The visibility function shows that the stellar contri- 
bution to the total 2.11 flux is less than ^ 40%, indicating a 
rather large optical depth of the circumstellar dust shell. The 2- 
dimensional Gaussian visibility fit yields a diameter of AFGL 
2290 at 2 . 1 1 fim of 43 mas x 5 1 mas, which corresponds to a di- 
ameter of 42 AUx50 AU for an adopted distance of 0.98 kpc. 

Our new observational results provide additional con- 
straints on the CDS of AFGL 2290, which supplement the in- 
formation from the spectral energy distribution (SED). To de- 
termine the structure and the properties of the CDS we have 
performed radiative transfer calculations for spherically sym- 
metric dust shell models. The observed SED approximately at 
phase 0.2 can be well reproduced at all wavelengths by a model 
with Toff = 2000 K, a dust temperature of 800 K at the inner 
boundary ri, an optical depth Ty ~ 100 and a radius for the 
single-sized grains of = 0.1 /im. However, the 2.11 fim 
visibility of the model does not match the observation. 

Exploring the parameter space, we found that grain size is 
the key parameter in achieving a fit of the observed visibility 
while retaining the match of the SED, at least partially. Both the 
slope and the curvature of the visibility strongly constrain the 
possible grain radii. On the other hand, the SED at longer wave- 
lengths, the silicate feature in particular, determines the dust 
mass loss rate and, thereby, restricts the possible optical depths 
of the model. With a larger grain size of 0.16 /^m and a higher 
TV = 150, the observed visibility can be reproduced preserving 
the match of the SED at longer wavelengths. Nevertheless, the 
model shows a deficiency of flux at short wavelengths, which 
is attributed to the model assumption of a spherically symmet- 
ric dust distribution, whereas the actual structure of the CDS 
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around AFGL 2290 is in fact non-spherical. Our study demon- 
strates the possible limitations of dust shell models which are 
constrained solely by the spectral energy distribution, and em- 
phasizes the importance of high spatial resolution observations 
for the determination of the structure and the properties of cir- 
cumstellar dust shells around evolved stars. 
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1. Introduction 

AFGL 2290 (OH 39.7H-1.5, IRAS 18560H-0638, VI 366 Aql) 
belongs to the group of type II OH/IR stars, which can be de- 
fined as infrared point sources with a maximum of the spectral 
energy distribution (SED) around 6 — 10 //m, with the 9.7 
silicate band in absorption, and with OH maser emission in the 
1612 MHz line (Habing |1996| ). Most of these objects show a 
long-period variability in the infrared and the OH maser emis- 
sion (Engels 1982; Herman & Habing 1985), although also a 



small fraction either varies irregularly with small amplitude or 
does not vary at all. OH/IR stars are surrounded by massive 
circumstellar envelopes composed of gas and small solid par- 
ticles (dust, grains). These circumstellar dust shells (CDS) are 
produced by the ejection of matter at large rates {M > 10~^ 
Mq yr^^) and low velocities (^ 15 kms^^), and in some cases 
they totally obscur the underlying star. Based on the luminosi- 
ties (-- 10"^ Lq), the periods (500d to 3000d) and bolometric 
amplitudes 1 mag), the kinematical properties and galac- 
tic distribution, the majority of OH/IR stars are highly evolved 
low- and intermediate-mass stars populating the asymptotic 
giant branch (AGB) (Habing 1996 1. They extend the sequence 
of optical Mira variables to longer periods, large r optical depths 
and higher mass l oss rates (Engels et al 
Lepine et al. |1995|). 



1983; Habing 1990 



The improvements of the observational techniques, espe- 
cially at infrared wavelengths, and the elaboration of increas- 
ingly sophisticated theoretical models have provided a wealth 
of new information on the structure, the dynamics, and the evo- 
lution of the atmospheres and circumstellar shells of AGB stars. 
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although many details still remain to be clarified (see the review 
by Habing 1996). A general picture has become widely ac- 



cepted in which both the large amplitude pulsations and the ac- 
celeration by radiation pressure on dust contribute to the mass 
loss phenomenon for AGB stars. From observations, correla- 
tions are found between the period and the infrared excess (in- 
dicating the mass loss rate) (DeGioia-Eastwood et al. 1981 



Jura 1986), and between the period and the terminal outflow 
velocity (Heske 1990| ). On the theoretical side, hydrodynami- 
cal models showed that due to the passage of shocks generated 
by the stellar pulsation the atmosphere is highly extended, thus 
enabling dust formation and the su bsequent acceleration of the 
matter (Wood 1979; Bowen 1988). The inclusion of a detailed 



treatment of dust formation revealed a complex interaction be- 
tween pulsation and dust formation, which results e.g. in a lay- 
ered dust distrib ution and af fects the derive d optical app earance 
(Fleischer et al. |1992| , |1995| ; Winters et al. |1994 |1995| ). 

Until now most interpretations of observations as well as 
most theoretical models are based on the assumption of a spher- 
ically symmetric dust shell, often motivated by the circularity 
of the OH maser maps. However, observations show that some 
objects have substantial deviations from spherical symmetry 
(e.g. Dyck et al. |1984]; Kastner & Weintraub |1994|; Weigelt 



et al. 1998). This suggests that the asymmetries obser ved in 



many post-AGB objects and planetary nebulae (cf. Iben |1995| ) 
may already start to develop during the preceding AGB phase, 
which provides new challenges for the modeling of the mecha- 
nisms and processes determining the structure of the dust shells 
around AGB stars. 

High spatial resolution observations can yield direct infor- 
mation on important properties of the dust shells around AGB 
stars, such as the dimensions and geometry of the shell. There- 
fore, such observations contribute additional strong constraints 
for the modeling of these circumstellar environments, which 
supplement the information from the spectral energy distribu- 
tion. Measurements of the visibility at near-IR wavelengths, for 
example, can be used to determine the radius of the onset of 
dust formation as well as to constrain the dominant grain size 



(Groenewegen 1997 1. To gain information on details of the spa- 



tial structure, in particular on asymmetries and inhomogeneties 
of the CDS, the interferometric imaging with large single-dish 
telescopes is especially well suited because one observation 
provides all spatial frequencies up to the diffraction limit of the 
telescope and for all position angles simultaneously, allowing 
the reconstruction of true images of the object. 

We have chosen AFGL 2290 for our study because it rep- 
resents a typical obscured OH/IR star with a high mass loss 
rate, whose location is not too far away from us. The distance 
to AFGL 2290 can be determined directly with the phase lag 



method (cf. Jewell et al. 1979), which gives D — 0.98 kpc (van 



Langevelde et al. 1990). For the bolometric flux at earth a value 
of /b - 2.4 IQ-i" Wm~2 jg derived by van der Veen & Rugers 



( 1989 ) from infrared photometry between 1 /im and 12 /im and 
the IRAS fluxes. At 0.98 kpc the luminosity is L = 7200 L©, 
which is within the typical range for an oxygen-rich AGB star 
The long period of P = 1424 d determined from the variation 



of the OH maser (Herman & Habing 1985) and the high mass 
loss rate suggest, that the star is in a late phase of its AGB evo- 
lution. 



So far. Chapman & Wolstencroft (1987) reported the only 
high angular-resolution infrared observations of AFGL 2290. 
From 1-dimensional slit-scan speckle interferometry with the 
UKIRT 3.8 m telescope at 3.8 /im and 4.8 /im they derive 1- 
dimensional visibilities and determine Gaussian FWHM di- 
ameters. Radiative transfer models for the AFGL 2290 dust 



shell have been presented by Rowan-Robinson ( 1982), Bedijn 



(1987), Suh (1991) and recently by Bressan et al. (1998). These 



models yield dust shell properties within the typical range of 
OH/IR stars, e.g. a dust mass loss rate of about 4 10~^ M,^ -"-"^ 
Bressan et al. 



(Bedijn 1987 



9.7 of about 10 (Bedijn 1987; Suh 1991). However, none 



1998), or an optical depth at 



of these studies includes constraints from high spatial resolu- 
tion infrared measurements. 

In Sect. 1^ we present the results of our speckle masking 
observations of AFGL 2290. The approach for the radiative 
transfer modeling is described in Sect. ^ comprising a short 
description of the code, the selection of the photometric data 
and a discussion of input parameters for the models. In Sect. Q 
we present the results of the radiative transfer modeling starting 
with the discussion of a model, which yields a good fit of the 
observed SED at all wavelengths but does not reproduce the ob- 
served 2.11 visibility. In search of an improved model the 
changes of the resulting SED and visibility under variations of 
the input parameters are investigated in the following sections. 
We finish the paper with a summary of the results and our con- 
clusions in Sect. a. 



2. Speckle masking observations 

The AFGL 2290 speckle data presented here were obtained 
with the Russian 6 m telescope at the Special Astrophysical 
Observatory (SAO) on June 14 and 16, 1998. We recorded a 
total number of 1200 speckle interferograms of AFGL 2290 
(600 on June 14 and 600 on June 16) and 2400 speckle inter- 
ferograms of the unresolved reference star HIP 93260 (1200 
on each of the two nights) with our 256x256 pixel NIC- 
MOS 3 camera through an interference filter with center wave- 
length 2.11 /im and a bandwidth of 0.192/^m. The exposure 
time per frame was 100 ms, the pixel size was 30.61 mas and 
the field of view 7'.'8 x 7'.'8. The 2.11 /im seeing was about 
^ 1'.'2. A diffraction-limited image of AFGL 2290 was re- 
constructed from the speckle data by the speckle masking bis- 



pectrum method (Weigelt 1977; Lohmann et al. 1983; Weigelt 



1991). The process includes the calculation of the average 



power spectrum and of the average bi-spectrum and the sub- 
traction of the detector noise terms from those. The modulus of 
the object Fourier transform was determined with the speckle 
interferometry method (Labeyrie 1970). The Fourier phase was 
derived from the bias-compensated average bispectrum. 

Figures |l] and ^ show the visibility function of AFGL 
2290 at 2.11 /im. The azimuthally averaged visibility decreases 
steadily to values below ^ 0.40 of the peak visibility at the 
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Fig. 1. Two-dimensional 2.11 /im visibility function of AFGL 
2290 derived from the speckle interferograms. The contour lev- 
els ai-e plotted from 20% to 80% of the peak value in steps of 
10%. 



Fig. 3. Diffraction-limited 2.11 /im speckle masking image of 
the AFGL 2290. North is at the top and east to the left. The 
contours level intervals are 0.25 mag. The lowest contour level 
is 3.25 mag fainter than the peak intensity. 




4 6 8 10 

Spatial frequency q [arcsec"^] 

Fig. 2. Azimuthally averaged 2.11 /zm visibility function of 
AFGL 2290 and errorbars. 



E 0.4 




40 60 80 

Distance from center [mas] 



Fig. 4. Azimuthally averaged image of AFGL 2290 (solid line) 
and of the unresolved reference star HIP 93260 (dashed hne). 



diffraction cut-off frequency (13.5 arcsec"^). Thus, the cir- 
cumstellar dust shell is almost totally resolved, and the con- 
tribution of the unresolved stellar component to the monochro- 
matic flux at 2.11 /im must be less than 40%, suggesting 
a rather high optical depth at this wavelength. In order to de- 
rive diameters for the dust shell, the object visibility function 
was fitted with an elliptical Gaussian model visibility function 
within a range of 1.5 arcsec"^ up to 7.5 arcsec^^. We obtain a 
Gaussian fit diameter of 43 mas x 5 1 mas for AFGL 2290 corre- 
sponding to 42 AUx50 AU for an adopted distance of 0.98 kpc 
or 5.7 x 6.8 r» for an adopted distance of 0.98 kpc and an 



adopted stellar radius of = 7.5 mas (cf. Sect. 4.1 ), respec- 
tively. 

Figure |^ shows the reconstructed 2.11 /im speckle mask- 
ing image of AFGL 2290. The resolution is 75 mas. Figure ^ 
shows the azimuthally averaged images derived from the re- 
constructed 2-dimensional images of AFGL 2290 and the ref- 
erence star HIP 93260. In the 2-dimensional AFGL 2290 image 
a deviation from spherical symmetry can be recognized. The 
intensity contours are elongated in the south-eastern direction 
along an axis with a position angle of 130°. 
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3. The radiative transfer modeling approach 

3.1. The radiative transfer code 

The radiative transfer calculations are performed with the code 



DUSTY developed by Ivezic et al. (1997), which is publicly 



available. The program solves the radiative transfer problem 
for a spherically symmetric dust distribution around a central 
source of radiation and takes full advantage of the scaling prop- 
erties inherent in the formulation of the problem. The formu- 
lation of the radiative transfer problem, the model assumptions 
and the scaling properties are described in detail by Ivezic & 
Elitzur ( 1997 ). Therefore, we give only a brief discussion here. 
The problem under consideration is a spherically symmetric 
dust envelope with a dust free inner cavity surrounding a cen- 
tral source of radiation. This geometry is not restricted to the 
dust shell of a single star. It can as well describe a dust enve- 
lope around a group of stars (e.g a binary) or even around a 
galactic nucleus. The radial dependence of the dust density be- 
tween the inner and outer boundary can be chosen arbitrarily. 
To arrive at a scale invariant formulation two assumptions are 
introduced: i) the grains are in radiative equilibrium with the 
radiation field, and ii) the location of the inner boundary ri of 
the dust envelope is controlled by a fixed temperature Ti of the 
grains at r i . Due to radiative equilibrium this temperature is de- 
termined by the energy flux at ri, which in turn is controlled by 
the energy flux from the central source via the radiative transfer 
through the dusty envelope. Then prescribing the dust temper- 
ature at 7'i is equivalent to specifying the bolometric flux at the 
inner boundary, and the only relevant propert y of t he input ra- 
diation is its spectral shape (Ivezic & Elitzur 1997 ). Similarly, 
if the overall optical depth of the dust envelope at some refer- 
ence wavelength is prescribed, only dimensionless, normalized 
distributions describing the spatial variation of the dust density 
and the wavelength dependence of the grain optical properties 
enter into the problem. 

This formulation of the radiative transfer problem for a 
dusty envelope is well suited for model fits of IR observations, 
because it minimizes the number of independent model param- 
eters. The input consists of: 

- the spectral shape of the central source of radiation, i.e. the 
variation of the normaUzed monochromatic flux with wave- 
length, 

- the absorption and scattering efficiencies of the grains, 

- the normalized density distribution of the dust, 

- the radius of the outer boundary in units of the inner bound- 
ary, 

- the dust temperature at the inner boundary, 

- the overall optical depth at a reference wavelength. 

For a given set of parameters, DUSTY iteratively determines the 
radiation field and the dust temperature distribution by solving 
an integral equation for the energy density, which is derived 
from a formal integration of the radiative transfer equation. 
For a prescribed radial grid the numerical integrations of radial 
functions are transformed into multiplications with a matrix of 
weight factors determined purely by the geometry. Then, the 



energy density at every point is determined by matrix inver- 
sion, which avoids iterations over the energy density itself and 
allows a direct solution of the pure scattering problem. Typi- 
cally fewer than 30 grid points are needed to achieve a relative 
error of flux conservation of less than 1 %. The number of points 
used in angular integrations is 2-3 times the number of radial 
grid points, and the build-in wavelength grid has 98 points in 
the range from 0.01 fim to 3.6 cm (see Appendix C in Ivezic & 



EHtzur 1997) 



The distributed version of the code provides a variety of 
quantities of interest including the monochromatic fluxes and 
the spatial intensity distribution at wavelengths selected by the 
user, but not the corresponding visibilities. Since we want to 
employ the visibilities obtained from our high spatial resolution 
measurements as constraints for the radiative transfer models, 
we have supplemented the code with routines for the calcula- 
tion of synthetic visibility functions. 

3.2. Selection of photometric data 

An important ingredient for the radiative transfer modeling of 
circumstellar dust shells around evolved stars is the spectral 
energy distribution (SED). Due to the variability of Miras and 
OH/IR stars, the SED of such objects ideally has to be deter- 
mined from coeval observations covering all wavelengths of 
interest. Unfortunately, no such coeval photometric data set for 
the wavelength region from A « 1 /im to A > 20 /im is avail- 
able in the literature for AFGL 2290. Thus, we have to define 
a 'composite' SED, which is derived from observations made 
by different authors at different epochs, but at about the same 
photometric phase (Griffin 1993 ). 

From the infrared photometry of AFGL 2290 available in 
the literature, we selected those publications which specify the 
date of observation and present the fluxes in tabulated form, 
either in physical units (e.g. Jy) or in magnitudes with given 
conversion factors (at least as a reference). Table |l] lists the ref- 
erences, the date and phase of observation and the wavelengths. 
The phases were determined from the period P = 1424d and 
the epoch of maximum, JD = 244 4860.8, which has been de- 
rived from the monitoring of the OH maser emission by Her- 
man & Habing (1985). Engels et al. (1983) determined periods 
of OH/IR stars from infrared observations and found that the 
periods and phases are in agreement for objects in common 
with the sample Herman & Habing (1985). 

It can be seen from the entries in Table [l| that the wave- 
length range from A = 1.2 /im to A = 20 /im is only fully 
covered by observations around phase 0.2 (see entries 2, 4, 9, 
12). The respective fluxes are shown in Fig. 5| 

The mea surements of Herman et al. (1984) and Nyman 
et al. (1993) match each other quite well at A = 3.8 and 
A = 4.8 /im, although the observations are separated by two 



peri ods. T he fluxes of Price & Murdock (1976) and Gehrz et 
al. (1985) agree with the Herman et al. and Nyman et al. data 
within the errors given by the authors. 

To represen t the S ED of AFGL 2290 we adopt the data of 
Herman et al. ( |l984| ), Gehrz et al. ( |l985| ), and Nyman et al. 
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Julian Date 


Phase* 


Ref. 










Wavelen 
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244 0000+ 


P = 1424 d 








1 


1045 


0.320 (-3) 


1 








4.2 








11.0 




19.8 


2 


2295 


0.198 (-2) 


1 
















11.0 




19.8 


3 


2725 


0.500 (-2) 


2 




2.2 


3.6 


5.0 


8.4, 8.8 




10.4 10.6 


11.6 


12.6 




4 


3726 


0.203 (-1) 


3 




2.3 


3.6 


4.9 


8.7 




10.0 


11.4 


12.6 


19.5 


5 


3972 


0.376 (-1) 


4 


1.25, 1.65 


2.2 


3.7 


4.8 














6 


4082 


0.453 (-1) 


3 




2.3 


3.6 


4.9 


8.7 




10.0 


11.4 


12.6 


19.5 


7 


4348 


0.640 (-1) 


3 




2.3 


3.6 


4.9 


8.7 




10.0 


11.4 


12.6 


19.5 


8 


4533 


0.770 (-1) 


4 


1.25, 1.65 


2.2 


3.7 


4.8 


8.2 


9.6 


10.2 




12.2 


19.6 


9 


5146 


0.200 (+0) 


5 






3.8 


4.8 


8.7 


9.7 


10.5 


11.5 


12.5 


20 


10 


7816 


0.075 (+2) 


6 


1.63 


2.23 


3.79 
















11 


7832 


0.087 (+2) 


7 


1.26, 1.68 


2.28 


3.80 
















12 


8041 


0.233 (+2) 


8 


1.24, 1.63 


2.19 


3.79 


4.64 















References: 1) Price & Murdock 



197f, 2) Lebofsky et alJ197f, 3) Gehrz et al. 1985, 4) Engels 1982, 5) Herman et al 



6) Noguchi et al. 1992 , 7) Xiong et al. 1994 , 8) Nyman et al 



1993 



* Numbers in parantheses give the cycle with respect to epoch JD 244 4860.8. 
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Fig. 5. IR-fluxes of AFGL 2290 observed about phase 0. 2 and 
at phase 0.77. The data are taken from Price & Murdock 1976 
(+), Her man et al. |1984| (o), Gehrz et al. |1985| (A), and Ny- 
man et al. 1993 (□). Also shown are the colour corrected IRAS 



fluxes adopted from van der Veen et al. 1995 (filled sqares), and 
the IRAS low resolution spectra (dashed line). The IRAS data 
are multiplied by a factor of 1 .95 in order to match the photo- 
metric data at 12 /im. The insert shows mm measurements by 



Walmsley et al. 1991 (•) and van der Veen et al. 1995 [SJ, 3 cr 
upper limits). 



( 1993 ). The scatter between the different data sets gives a rough 
estimate of the uncertainty of the 'composite' SED at phase 0.2 
of w 0.25. We do not correct for interstellar extinction because 
the corrections are less than, or of the same order as, the uncer- 



tainty estimated above. For AFGL 2290 Herman et al. (1984) 
give a value of Ay = 1.6 for an adopted distance of 1.19 kpc 
which reduces to Ay ~ 1.3 at a distance of 0.98 kpc. With 
the wavelength dependence of the interstellar extinction in the 



infrared from Becklin et al. ( 1978 1 one obtains a correction fac- 
tor of 1.35 at A = 1.25 /im, 1.03 at A = 4.8 /im, and 1.15 at 
A = 9.5 /im. 

AF GL 22 90 was observed by IRAS (IRAS Point Source 
Catalog 1985). We adopt the colour corrected broadband fluxes 
given by van der Veen et al. (1995) and the IRAS low resolu- 
tion spectra from the IRAS Catalog of Low Resolution Spec- 
tra (1987). The latter are corrected according to Cohen et al. 
(1992). Since the broadband fluxes and spectra are averages 



of several measurements taken at different phases (IRAS Ex- 



planatory Supplement 1985), the flux levels, e.g. at 12 /im, are 
lower than the fluxes from ground based observations around 
phase 0.2. Therefore, we multiply the IRAS data with a factor 
of 1 .95 to join them with the ground based data. 

Finally, observations at mm wavelengths were reported by 
Walmsley et al. ( 1991 ) who measured a flux of 0.025 Jy at A = 
1.25mm at phase 0.18 (JD 2447960), and by van der Veen et 
al. ( 1995 ) who derived 3(t upper limits of 0. 1 3 Jy and 0. 14 Jy at 
0.76 mm and 1 . 1 mm respectively for phase 0.25 (JD 2448069), 
which are consistent with the 1.25 mm flux. 

3.3. Selection of input parameters 

We represent the central star by a blackbody with an effective 
temperature Teff- In contrast to the visible M type Mira vari- 
ables with Toff ^ 3500 K, the effective temperature of OH/IR 
stars with optically thick dust shells cannot be directly deter- 
mined. However, if OH/IR stars can be considered as an ex- 
tension of the Mira sequence to longer periods and larger opti- 
cal depths, one might extrapolate the period-Tcff relation for 



Mira variables derived by Alvarez & Mennessier (1997) to 
P > 650 d, which yields Tcff < 2500 K in agreement with 
the values expected for the tip of the AGB. 
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The dust density distribution is obtained from the velocity 
law, which results from an approximate analytic solution for a 
stationary dust driven wind with constant mass loss rate (e.g. 



Schutte & Tielens 1989). If the gas pressure force is neglected 



and the flux averaged absorption coefficient is assumed to be 
constant with the radius r in the wind, the velocity distribution 
is given by 



v{r) 



1 - 




(1) 



where vi denotes the velocity at the inner boundary ri, and Voo 
is the velocity at infinity. The relevant free parameter is the ratio 
of these velocities S = vi/voo, because only the normalized 
density distribution enters into the calculation. 

We adopt this velocity law, because it accounts for the 
changing density gradient due to the acceleration of the matter 
by radiation pressure on dust in the innermost parts of the dust 
shell. Compared to a dust shell with al/r^ density distribution 
and equal optical depth the dust density at ri is higher by a fac- 
tor of 0.5(1 + (5) /(5 and the mass loss rate is lower by a factor of 
0.5(1 + (5) (cf. Le Sidaner & Le Bertre |1993[ ). According to the 
theory of dust driven winds the velocity at the inner boundary, 
where efficient grain condensation takes place and the acceler- 
ation of the matter by radiation pressure on dust starts, is close 
to local sound velocity Cg (see Gail 1990). This is supported 



by observations of the velocity separation of the SiO maser 



emission in OH/IR stars (Jewell et al. 1984 ), which presumably 
originates from the dust forming region. With Cg < 2kms^^ 
for temperatures of about 1000 K and with the measured out- 
flow velocity of AFGL 2290 of v^o — 16 kms^^ (Herman & 
Habing 1985) one obtains S w 0.12, which we adopt as the 
standard value for 5. 

As described in the previous section, the location of the in- 
ner boundary ri of the dust shell is determined by the choice of 
the dust temperature Ti at ri . For the outer boundary rout we 
adopt a default value of rout = 10"^ ri. As shown in the follow- 
ing section a larger outer boundary affects only the far infrared 
fluxes for A = 100 fim without altering the other properties of 
the model. 

We consider spherical grains of equal size described by 
the grain radius Ogi . This is certainly a simplification because 
based on theoretical and observational arguments, one expects 
the presence of a grain size distribution n(agi ). Therefore, 
a size distribution similiar to the one observed in the ISM 
(n(agr) oc a~^^^) is often assumed for radiative transfer mod- 



els of circumstellar dust shells (e.g. Justtanont & Tielens |1992| , 
Griffin 



(Gauger et al. 199Cl| , Winters et al. 1997) which might roughly 
be approximated by a single dominant grain size. 

For the dust optical properties we adopt the complex 
refractive index given by Ossenkopf et al. (1992) for 'warm, 
oxygen-deficient' silicates. The authors consider observational 
determinations of opacities of circumstellar silicates as well 
as laboratory data and discuss quantitatively the effects of in- 
clusions on the complex refractive index, especially at shorter 
wavelengths (A < 8 /im). These constants yield a good match 
of the overall spectral shape of the observed SED of AFGL 
2290, especially of the 9.7 /im silicate feature. However, we 
will also discuss the effects on the radiative transfer models 
resulting from different optical constants in the Appendix. 
With the tabulated values of the complex refractive index the 
extinction and scattering efficiencies are calculated from Mie 
theory for spherical particles assuming isotropic scattering. 

Once a satisfactory fit of the spectral shape is achieved with 
suitably chosen values for the remaining input parameters Tcs, 
Ti, Ogr, and ro.55 (the optical depth at the reference wavelength 
0.55 /Ltm) the match of the normalized synthetic SED with the 
observed SED determines the bolometric flux at earth /b. Com- 
bined with the effective temperature one obtains the angular 
stellar diameter 9^, and thereby the spatial scale of the system. 
With an assumed distance D to the object, the luminosity L«, 
the radius of the inner boundary in cm, and the dust mass loss 
rate Md can be calculated. The latter quantity is given by: 



Md = 27rriWoo(l + <5)pgr-~ 



Td(A) 



t(A) 



(2) 



where denotes the specific density of the grain material, 
Td(A) is the dust optical depth at wavelength A, and Qcxt = 
Kcxt/V"gr is the extinction cross section Koxt per unit volume 
Vgr of the grains. For single sized grains Qoxt is proportional to 
the extinction efficiency divided by the grain radius Qcxt /«gr> 
and in the Rayleigh limit 27ragr ^ A it is independent of the 
grain radius Ogr. For the specific density of silicate grains we 
adopt Pgr = 3gcm~'^ as a typical value. For the distance to 
AFGL 2290 we use D = 0.98 kpc (van Langevelde et al. 
19901). 



4. Radiative transfer modeling of AFGL 2290 

4.1. A radiative transfer model for the SED of AFGL 2290 
Starting with the parameters of previous radiative transfer mod- 



els for AFGL 2290 presented by Rowan-Robinson (1982), Be 



1993). Consistent models for stationary dust driven dijn (1987 1, and Suh (1991 ), we achieved a satisfactory match 



winds, which include a detailed treatment of (carbon) grain 
formation and growth, in fact yield a broad size distribution 
which can well be approximated by a power law (Dominik 
et al. 1989| ). However, in circumstellar shells around pulsat- 
ing AGB stars the conditions determining the condensation of 
grains change periodically. The time available for the growth of 
the particles is restricted by the periodic variations of the tem- 
perature and density. This results in a narrower size distribution 



of the observed SED with the set of parameters given in Table 
^ after a few trials. Henceforth we will refer to this parameter 
set as model A, and we will first discuss its properties before 
we use it as a reference for an investigation of the sensitivity of 
the results on the parameters. 

The SED of model A is shown in Fig. || Figure ^ displays 
an enlargement of the 5 — 2 5 /im region with the 9 . 7 /im silicate 
feature. From the shortest wavelength at A = 1.25 /im up to 
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Table 2. Parameters and resulting properties of model A. 



Teff [K] 


Ti [K] 


Ugr [/im] 


Tout/ri 


To. 55 


2000 


800 


0.10 


10* 


100 


Md [Meyr-i] 


ri [R*] 


/b [Wm-2] 


[mas] 


TlO 


2.7 10"'' 


7.80 


3.0 10"^° 


7.50 


7.49 



A = 1.25 mm model A provides a good fit to the observations. 
The location, shape and strength of the silicate feature around 
10 /im is well reproduced with the adopted optical data from 
Ossenkopf et al. ( 1992 ). Only in the 18 /im region is there a 
noticeable deviation because the model shows a weak, broad 
emission which is absent in the IRAS LRS spectrum. 

In addition to the input parameters. Table || also lists the 
derived properties of model A. Our value for the bolomet- 

^ is consistent 



-10 Wm- 

"^"Wm"^ determined by van der 

Veen & Rugers (|l989|). With Toff = 2000 K the angular stel- 



ric luminosity at earth of /b — 3 10 
with the values of 2.4 10" 



lar diameter is 9^, = 7.5 mas and we obtain a stellar radius 
i?* = 790 R0 and a luminosity = 9000 L© adopting a dis- 
tance D — 0.98 kpc. The bolometric flux should be quite ac- 
curate considering the quality of the fit. However, the errors in 
the distance determination (e.g. an — 0.34, van Langevelde et 
al. 1990| ) and the uncertainty in the determination of Tcs from 
the radiative transfer modeling (see Sect. 4.2.2 ) are rather large, 
resulting in correspondingly large uncertainties for L*, i?* and 

The derived dust mass loss rate of 2.7 10"'' M0yr^^ is 
close to the values derived fr om ra diative transfer models by 
other authors. Bressan et al. (|1998| ) obtain 4.5 lO"'^ Mqjt-^, 
and the model of Bedijn ( |1987| ) yields|^ 4 lO"'^ Mgyr^^. From 
a relation between the strength of the 10 /^m feature and the 
color temperature Schutte & Tielens ( |1989| ) obtain ed — 



2.4 10"^ Mgyr-i for AFGL 2290, and Heske et al. ( |1990t ) es- 



timated Md = 1.2 10"'^ Mgyr^i from the 60 /im IRAS flux. 

All together the parameters and derived properties of model 
A lie in the typical range of values obtained from radiative 
transfer models for OH/IR stars showing the silicate feature in 
absorp tion. In the calculations of Lorenz-Martins & de Araujo 
(|1997|) Toff ranges from 1800 K to 2400 K, Ti from 650 K to 



1200 K, and T9.7 from 7 to 17. Justtanont & Tielens ( |1992D 
derive dust mass loss rates between 2.6 10"^ Moyr~^ and 
2.21O~''M0yr~\ and optical depths T9.7 between 4.85 and 
19.6 for their sample of OH/IR stars. 

In addition to model A with rout It\ = 10*, Fig. || also dis- 
plays the SED for a model with the same parameters but with an 
outer boundary, which is ten times smaller, i.e. rout 1^1 — 10"^ 
(dotted line). The spectra are virtually indistinguishable up to 
A ^ 60/im. The additional cold dust due to the larger outer 
boundary of model A increases the far infrared fluxes (see in- 
sert of Fig. H), but does not affect the SED at shorter wave- 



lengths (cf. Justtanont & Tielens 1992). Since the values of the 



' We have calculated from the values for ri and rio given by 
Bedijn \\9i% adopting Qcxt(10 /im) = 5.6 10^ cm"\ 
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10 

Wavelength [urn] 



100 



Fig. 6. Spectral energy distribution for model A (solid line) and 
adopted photometry of AFGL 2290 at about phase 0.2 (see Fig. 
H for the references and corresponding symbols). The insert 
shows the SED in the mm range. The dotted line displays the 
SED for a model with the parameters of model A and an outer 
boundary of rout = lO'^ ri. 



„-r- 1-5 



E 




15 

Wavelength [|j,m] 

Fig. 7. 5 — 25 /im spectrum of model A (solid line) and the 
corrected IRAS LRS spectra (dashed lines). The IRAS data are 
scaled by a factor of 1.95 to match the ground based photome- 
try (see Fig. S for the references and corresponding symbols). 



derived properties (/b, M^, ■■■) are determined essentially by 
the shape of the SED below 60 /im, they do not change for 
Tout/ri ^10^ 

Figure II shows a comparison of the 2.11 /im visibility cal- 
culated for model A with the azimuthally averaged 2.11 /im 
visibility Vohs derived from our speckle observations. Although 
model A yields a good fit to the observed SED, it fails to repro- 
duce the observed visibility. The slope of the model visibility 
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Table 3. Resulting properties for models with the parameters 
of model A, but with different dust temperatures at the inner 
boundary. 





To. 55 


Md 






no 


,/b 




[K] 




[Moyr- 




[R*] 




[Wm- 




600 


65 


2.86 10- 


-7 


12.89 


4.87 


2.60 10 


-iO 


800 


100 


2.66 10" 


-7 


7.80 


7.49 


3.00 10 


-10 


1000 


150 


2.70 10" 


-7 


5.27 


11.23 


3.12 10 


-10 


1200 


200 


2.66 10" 


-7 


3.91 


14.97 


3.25 10 


-10 



4 6 8 10 

Spatial frequency q [arcsec"^] 



12 



Fig. 8. Synthetic 2.11 /im visibility for model A (solid line) 
and the azimuthally averaged 2.11 /im visibility obtained from 
the speckle observation (squares). 



is much shallower than the slope of the observed one, and the 
model visibility levels off at g ^ 20 arcsec^^ with V2.11 = 0.6. 
From the observed visibility we obtained an upper limit for the 
contribution of the star to the monochromatic flux at 2.11 /im 
of « 40 %, whereas model A yields 62 %, indicating that the 
model optical depth of T2.11 = 3 at 2.11 /im is too low. Fur- 
thermore the steeper slope of Vohs suggests that a levelling- 
off should occur at a smaller spatial frequency compared to the 
model visibility, which means that the observed intensity distri- 
bution would be rather more extended than the intensity profile 
of model A. 

To summarize, a dust shell model, which yields a visibility 
in agreement with the observations, requires a larger optical 
depth at 2.11 /im and a more extended intensity distribution 
compared to model A, but it has to produce the same spectral 
energy distribution. It turns out that these requirements can be 
partially fulfilled with suitably modified parameter values. 

4.2. Effects of parameter variations on the calculated SED 
and the model visibility 

In order to find a model which matches both the observed SED 
and the visibility, as well as exploring the sensitivity of the 
SED and the 2.11 /xm visibility of the models on the input 
parameters, we have calculated various sequences of models. 
In each sequence one parameter was varied within a certain 
range, while the other parameters were kept fixed at the val- 
ues of model A, except for the optical depth at the reference 
wavelength ro.55. This quantity was adjusted for each model in 
order to obtain a match of the SED with the observation, es- 
pecially with the observed strength of the 9.7 /im silicate fea- 
ture. It turns out that models, for which the SED fit the feature 
equally well have (almost) equal dust mass loss rates, i.e. with 
our procedure we force the models in a sequence to have equal 
mass loss rates instead of equal ro.55. 



4.2. 1 . Effects of different dust temperatures at the inner 
boundary 

The effects of different dust temperatures at the inner boundary 
on the calculated SED and the 2.11 /im visibility are displayed 
in Fig. ^. The dust temperature is varied between Ti — 600 K 
and Ti — 1200 K, and the derived properties of the correspond- 
ing models are given in Table |[ 

For the values of Ti presented here, the SED is changed 
at wavelengths below A w 9 /im. Compared to model A with 
Ti = 800 K a smaller value of Ti results in a higher flux below 
A = 3 /im and a deficiency of flux at wavelengths between 
3 /im and 9 /im. Values of Ti > 800 K produce a deficiency of 
flux below A = 3 /im and an excess of flux between 3 /im and 
9 /im. As a consequence, the bolometric fluxes are different for 
these models. 

The changes of the SED for a variation of Ti with fixed 
other parameters, including T0.55, have been studied already by 
Ivezic & Elitzur (1997). For high ro.55 > 100 decreasing Ti 



lowers the flux at shorter wavelengths, increases the strength 
of the silicate feature and raises the flux at larger wavelengths, 
similiar to the changes produced by increasing ro.55. Since we 
have adjusted ro.55 for each value of Ti to reproduce the ob- 
served strength of the 9.7 /im silicate feature, the model SED 
are only changed at shorter wavelengths in a modified way. 
This behaviour can be understood from the competition of the 
effects caused by lowering Ti and simultaneously decreasing 
''0.55 • Up to A ^ 3 /im the increase of the monochromatic flux 
induced by the smaller optical depth more than compensates 
the decrease of the flux due to a lower Ti, but in the region 
3 /im ^ A ^ 8 /im the effect of lowering Ti dominates and the 
monochromatic flux decreases with Ti. 

The variation of Ti mainly affects the slope of the 2.11 /im 
visibility V2.11 via the different ro.55. The optical depth in- 
creases with Ti resulting in a broader intensity profile and a 
smaller stellar contribution to the monochromatic flux. This 
leads to a steeper decline of V2.11 (see Ivezic & Elitzur 1996). 



The curvature of the visibility only noticeably changes at 
higher values of the spatial frequency. This change of the slope 
of the visibility indicates the onset of a levelling off of V2.11, es- 
pecially for the model with Ti = 600 K. Since its inner bound- 
ary is located at 12. 9i?* or 97 mas, the visibility should ap- 
proach a constant value at q w 13arcsec~^. An increase of 
Ti above 1000 K yields only a marginally steeper slope of the 
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Fig. 9. SED (a) and visibilities (b) for models with the pa- 
rameters of model A, but with different dust temperatures at 
the inner boundary: Ti ~ 600 K (long dashed line), 800 K 
(sohd line), 1000 K (short dashed Hne), 1200 K (dotted Hne). 
The corresponding optical depths and derived model properties 
are given in Table ^ The spectra have been scaled with differ- 
ent /b to match the observations at A > 8 fim. The symbols 
denote the observations (see Figs. ^ and^). 
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Fig. 10. SED (a) and visibilities (b) for models with the pa- 
rameters of model A, but with different effective temperatures: 
Teff = 1600 K (dashed line), 2000 K (solid Hne), 2400 K (dot- 
ted line). The corresponding optical depths and derived prop- 
erties are given in Table The spectra have been scaled with 
different /b to match the observations at A > 8 /im. The sym- 
bols denote the observations (see Figs. ^ and^). 



visibility, although ro.55 increases and the stellar contribution 
to the 2.11 /im flux correspondingly decreases. However, the 
inner boundary is shifted simultaneously to smaller radii, re- 
sulting in similiar slopes of V2.11 below q ^ 13 arcsec^^. 

To summarize, the changes of the SED due to a variation of 
Ti can in principle be compensated by a corresponding varia- 
tion of To. 55, but the presence of the silicate feature constrains 
the choice of the optical depth to values which reproduce the 
strength of the observed feature, i.e. to combinations of Ti 
and ro.55 which yield a fixed value for oc tiTio. Increas- 
ing Ti yields a steeper slope of the 2.11 /im visibility. Above 
Ti = 1000 K, however, these changes are small and will not 
result in a model which matches the observed visibility. 



4.2.2. Effects of different effective temperatures on the SED 
and the visibility 

Figure |l^ displays the effects of different effective tempera- 
tures on the calculated SED and the 2.11 /im visibility. The 
effective temperature is varied between Tcs = 1600 K and 
Toff — 2400 K. The derived properties of the corresponding 
models are given in Table | 

The effects of the variation of Toff on the model SED are 
qualitatively similiar to the effects induced by a variation of Ti. 
Lowering Toff from 2400 K to 1600 K decreases the monochro- 
matic flux below A ^ 2.5 /im and increases the flux at longer 
wavelengths up to A = 8 /im. The main cause for these changes 
is the shift of the wavelength X„i, where the stellar (black- 
body) spectrum reaches its maximum, from A„j = 1.5 /im to 
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Table 4. Resulting properties for models with the parameters 
of model A, but with different effective temperatures. 





Tf).55 


Md 




r-L 


TlO 


9, 


/b 




[K] 




[Moyr- 




[R,] 




[mas] 


[Wm- 




1600 


110 


2.61 10^ 


-7 


4.45 


8.23 


11.7 


3.110" 


-w 


2000 


100 


2.66 10" 


-7 


7.80 


7.49 


7.5 


3.0 10" 


-10 


2400 


92 


2.73 10" 


-7 


12.60 


6.86 


5.2 


2.910" 


-10 



A„j = 2.3 /im. The monochromatic flux of the 'hotter' star is 
larger below a certain wavelength, which is, in addition, af- 
fected by the slightly different optical depths. However, for the 
large optical depths considered here, the effects of different 
Toff on the SED are small as already shown by e.g. Rowan- 
Robinson (1980). It is interesting to note that the changes of 



the SED due to a variation of Ti can be compensated for by a 
corresponding variation of Tcs, at least within a certain range 
of values. 

The variation of T^s has only a negligible effect on the 
2.11 /im visibility for the temperature range considered here. 
Since the changes of the optical depths are small, both the an- 
gular diameter of the inner boundary and the stellar contribu- 
tion to the 2.11 fim flux only moderately increase with Toff- 
Correspondingly, the visibility approaches a slightly higher 
constant value at a slightly smaller spatial frequency result- 
ing in the minor differences for V2.11 at spatial frequencies 
q < 13.5 arcsec"^. Thus, changing the Tcs within a reasonable 
range cannot produce a model, which matches the observed 
visibility. 
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4.2.3. Effects of different grain radii 

The effects of different grain radii on the calculated SED and 
the 2.11 /Ltm visibility are displayed in Fig. [Tl|. The grain radius 
is varied between Ogr — 0.04 /im and Ogr — 0.12 /im. The de- 
rived properties of the corresponding models are given in Table 
||. Figure |l^ shows the extinction coefficient per unit volume of 
the grains Kcxt/V^i obtained from the optical data for 'warm' 
silicates from Ossenkopf et al. (1992). 



The choice of the grain radius only affects the short wave- 
length tail of the SED below A ^ 3 fim, because at these 
wavelengths Koxt/V^r still depends on agr, but it becomes in- 
dependent of flgr at longer wavelengths (see Fig. |l2|). This be- 
haviour is caused by two factors: the contribution of scattering 
to extinction and the dependence of the absorption efficiency 
on the grain size. The scattering efficiency per unit volume 
of the grains, which is cx Ogj., steeply declines with increas- 
ing wavelength and can be neglected above a certain wave- 
length depending on agr. The absorption efficiency depends on 
the grain size only at short wavelengths and becomes indepen- 
dent of agr once the grains are sufficiently small compared to 
the wavelength. Therefore, the grain radius is constrained by 
the observed fluxes at the shortest wavelengths A ^ 2 /tm. 
In our case the photometry at A = 1.65 /tm excludes grain 



Fig. 11. SED (a) and visibilities (b) for models with the param- 
eters of model A, but with different grain radii ranging from 
agr = 0.04 /im to agr = 0.16 /im. The corresponding optical 
depths and derived model properties are given in Table | The 
symbols denote the observations (see Figs. ^and^. 



radii agr ^ 0.12 /im and agr ^ 0.06 /im and the photome- 
try at A = 1.25 /tm restricts the grain radii to values close to 
agr = 0.1 /im. However, the values of the absorption and scat- 
tering efficiency depend on the adopted optical data. A different 
data set can result in vastly different grain radii (see Appendix 
A). 

The variation of the grain radius has two effects on the 
2.11 /im visibility. First, the slope of V2.11 steepens with in- 
creasing values of agr, because models with larger grain radius 
require a higher optical depth at 2.11 /im in order to match the 
observed SED for A > 2 /im. With increasing optical depth 
the intensity distribution becomes broader and the stellar con- 
tribution to the monochromatic flux at 2.11 /im decreases. Cor- 
respondingly, the decline of visibility with spatial frequency 



becomes steeper (see Ivezic & Elitzur 199(: ). The second effect 
is the change of the curvature of V2.11, which is noticable at 
low spatial frequencies. The curvature changes its sign at about 
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Table 5. Resulting properties for models with the parameters 
of model A, but with different grain radii. 



agr 


To. 55 


Md 




r\ 


no 


T"2.2 


[nm] 




[Mo yr^ 




VRA 






0.04 


20.0 


2.60 10" 




6.93 


8.22 


2.86 


0.06 


37.5 


2.69 10" 


7 


7.15 


8.24 


2.96 


0.08 


65.6 


2.66 10" 


7 


7.41 


7.86 


2.98 


0.10 


100 


2.66 10" 


7 


7.80 


7.49 


3.07 


0.12 


140 


2.63 10" 


7 


8.22 


7.03 


3.20 


0.14 


170 


2.60 10- 


7 


8.70 


6.59 


3.40 


0.16 


150 


2.65 10" 


7 


9.24 


6.33 


3.76 



/b = 310"'" Wm" 



= 7.5 mas. 



Ogr < 0.1 /im. This behaviour reflects the changes of the spa- 
tial intensity distribution. At large radial offsets h from the star 
the intensity decreases approximately as oc 6"^, because 
the optical depth al ong t he line of sight at b becomes small 
(see Jura & Jacoby 1976 ). For smaller offsets 6, however, the 
decline of the intensity steepens and the slope of the visibil- 
ity changes accordingly, i.e. the curvature of V indicates the 
changing slope of the spatial intensity distribution. Thus, the 
visibility constrains the grain radii not only via its slope, but 
also via its curvature. 

The observed visibility T4bs declines in almost a straight 
line to values below ^ 0.40 at g = IS.Sarcsec""^ with only a 
slight curvature. V^obs is fairly well matched by the model with 
Cgi- = 0.16 /im, although the curvature of the model visibil- 
ity is a little too strong. Since the visibilities for models with 
Ogr = 0.15 /im and Ogr = 0.17 /im already fall outside the er- 
ror bars of the observation, at least for certain spatial frequen- 
cies, the grain radius is deter mined by Vobs with an uncertainty 
< 0.01 /im (cf. Groenewegen 1997 ). 

However, the SED for the model with Ogr = 0.16 /<m, 
shows a deficit of flux below A = 3 /im. This deficit cannot be 
removed by a change of Ti or Toff. Although lowering Tf in- 
creases the flux at A < 3 /im it also decreases the flux at longer 
wavelengths. Furthermore, the inner boundary of the dust shell 
is moved outwards which deforms the resulting visibility in a 
way that destroys the fit. Increasing Tcff yields a similiar be- 
haviour. 

The deficit of the short wavelength model flux could have 
several causes, but the clear evidence for a non-spherical dust 
distribution around AFGL 2290 from our speckle masking ob- 
servations favors the explanation that the deficiency of flux in 
the model is due to the assumption of a spherically symmetric 
circumstellar dust shell. 

A more general assumption would be that the CDS has an 
axisymmetric, 'disk-like' structure. Theoretical investigations 
show that the variation of the effective optical depth with the 
inclination of a disk-like dust distribution affects the shape 
of the SED up to far infrared wavelengths, as well as the 
monochromatic intensity distributions and the c orresp onding 
visibilities (e.g. Efstathiou & Row an-Robinson 1990 ; Colli- 
son & Fix 1991 i Lopez et al. 1995 ; Men'shchikov & Henning 
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agr = 0.12 (im 
a„r = 0.14 urn 
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Fig. 12. Extinction coefficient per unit volume of the grains 
Kcxt/V^r for different grain radii ranging from Ogr = 0.04 /im 
to Ogr = 0.16 /im calculated with the optical data from Os- 
senkopf et al. (1992) for 'warm' silicates. 



1997 ). If the disk-like dust distribution is viewed at an interme- 
diate inclination one expects more flux at visual wavelengths 
than in the case of a spherically symmetric dust distribution, 
due to scattered light escaping from the optically thinner polar 
region located either above or below the equatorial plane. In 
other words, we expect a deficiency of the model flux at short 
wavelengths if we model the SED of an aspherical dust distri- 
bution assuming a spherically symmetric dust distribution. 

5. Summary and conclusions 

We have presented the first diffraction-limited 2.11 /im speckle 
masking observations of the circumstellar dust shell around the 
highly obscured type 11 OH/IR star AFGL 2290. The resolu- 
tion achieved with the SAO 6 m telescope is 75 mas, which is 
sufficient to partially resolve the circumstellar dust shell at this 
wavelength. From a 2-dimensional Gaussian fit of the visibil- 
ity function the diameter was determined to be 43 mas x 5 1 mas, 
which corresponds to a diameter of 42 AU x 50 AU for a dis- 
tance of 0.98 kpc. The reconstructed image shows deviations 
from a spherical structure with an elongation at position angle 
130°. 

Our new high resolution spatial measurements provide ad- 
ditional strong constraints for radiative transfer models for the 
dust shell of AFGL 2290, supplementing the information pro- 
vided by the spectral energy distribution (SED). In order to 
investigate the structure and the properties of the circumstel- 
lar dust shell we have performed radiative transfer calculations 
assuming a spherically symmetric dust distribution. The spec- 
tral energy distribution at phase ~ 0.2 can be well fitted at 
all wavelengths by a model with an effective temperature of 
2000 K, a dust temperature at the inner boundary of 800 K, an 
optical depth at 0.55 /J,m of 100, and a radius for the single- 
sized grains of 0.1 /im, using t he opt ical constants for 'warm' 
silicates from Ossenkopf et al. (1992). From this fit we derived 
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e.g. a bolometric flux at earth of 3 10^^° Wm~-^, a radius of 
the inner boundary of the dust shell of ri = 7.8 i?*, and a 
dust mass loss rate of 2.7 10~^ MQyr~^, in agreement with the 
results of previous radiative transfer models for AFGL 2290. 
However, this model does not reproduce the observed 2.11 /im 
visibility function. 

We have, therefore, investigated the changes of the calcu- 
lated SED and the model visibility with the input parameters in 
search of an improved model. We found that the grain size is the 
key parameter in achieving a fit of the observed visibility, while 
retaining at least a partial match of the SED. Both the slope and 
the curvature of the visibility react sensitively to the assumed 
grain radii. With the assumption of single-sized grains we ob- 
tain an uncertainty of less then ±0.01 /xm for Cgi . Another re- 
sult was that the dust mass loss rate is well constrained by the 
shape of the SED at longer wavelengths and, especially, by the 
shape of the silicate absorption feature. For given optical con- 
stants the value of the dust mass loss rate, as derived from the 
match of the feature, is not very sensitive to changes of the in- 
put parameters. The uncertainty of is ^ 3 10~* M0yr~^. 
The effective temperature and the dust temperature at the in- 
ner boundary, however, are less well constrained. We roughly 
estimate a range of ±300 K for Toff and ±100 K for Tf. 

The shape of the observed visibility and the strength of the 
silicate feature constrain the possible grain radii and optical 
depths of the model. The observed visibility can be reproduced 
by a model with a larger grain size of 0.16 /im and a higher 
Tv = 150, preserving the match of the SED at longer wave- 
lengths. Nevertheless, the model shows a deficiency of flux at 
short wavelengths, which can be explained if the dust distribu- 
tion is not sperically symmetric. If the CDS of AFGL 2290 has 
in fact a disk-like structure, the radial optical depths vary be- 
tween the equatorial and polar direction. Due to the scattered 
radiation escaping from the optically thinner polar regions one 
expects more flux at shorter wavelengths than from a spheri- 
cally symmetric system with equal optical depth towards the 
star. 
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Appendix A: Effects of different optical constants 



Table A.l. Resulting properties for models with the parameters 
of model A, but with different dust optical properties. 



In Sect. 4.2.3 the deficiency of short wavelength flux for the 
model which matches the observed 2.11 /im visibility of AFGL 
2290 was explained by the assumption of a spherically sym- 
metric dust distribution. However, the deficiency could be due 
to a different cause. For example, the optical constants of 'as- 
tronomical' silicates at wavelengths < 7.5 are not well 
known (Ossenkopf et al. 1992| ) and it is possible that the op- 
tical properties of the grains around AFGL 2290 at short wave- 
lengths differ from our assumption. 

Therefore, we have investigated the effects of different dust 
optical properties on the SED and the visibility, as shown in 
Fig. A.l. The models have been calculated with the parame- 
ters of model A (except for the value of ro.55) using the opti- 



cal data from Ossenkopf et al. (]1992) for 'cold' silicates (sil 



Oc), from Draine & Lee (1984) (sil-DL) and from David & 



per unit 
and 



A.2 



A.l 



Pegourie ( |1995| ) (sil-DP). The extinction coefficient 
volume of the grains for Ogr = 0.1 iim is shown in Fig. 
the derived properties of the models are given in Table 

The differences of the Kcxt/V^gr resulting from the optical 
data sets are more or less directly translated into modifications 
of the SED, if the different values for T0.55 are taken into ac- 
count Compared to the 'warm' silicates of Ossenkopf et al. 
(1992) (sil-Ow), the extinction of sil-Oc grains is higher be- 
tween A « 1.3 /im and 8 /im, resulting in a lower monochro- 
matic flux of the corresponding model. For the sil-DL and sil- 
DP data the extinction is lower resulting in an excess of flux. 
Because the shape of the silicate features at around 10 /im and 
18 /im is similiar for the sil-Ow, sil-Oc, and sil-DP data, ex- 
cept for a slightly different ratio of the peak strengths, they 
yield comparably good fits to width and strength of the ob- 
served feature, which has its center at 10 /im. In contrast, the 
silicate feature from the sil-DL data peaks at 9.7 nm, and it is 
broader than the observed one. 

Because the value of Afd depends on the adopted op- 
tical properties of the grains (see Eq. we obtain 
for the different models ran ging from 1.6 10~^ Moyr~^ to 
2.71O"'^M0yr~i (see Table |A.1[). Nevertheless, if the dust 
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mass loss rate is derived from the match of the silicate absorp- 
tion feature, its value is not very sensitive to variations of the 
effective temperature, dust temperature at the inner boundary 
and the grain radius as long as the models are calculated with 
the same optical constants. 

The changes of the 2.11 /im visibilities are again caused by 
the different optical depths of the models at this wavelength. 
For other fixed parameters a higher optical depth produces a 
more extended brightness distribution and, thereby, a steeper 
decline of the visibility. The optical depth at 2.11 /im has 
similiar values for the sil-Ow and sil-Oc models and lower, 
but again similiar values for the sil-DL and sil-DP models. 
Hence, the decline of visibilities from the latter models is 
shallower. 



The optical properties from David & Pegourie ( 1995 ) yield 
a fit of the silicate feature, which i s com parable to the fit ob- 
tained with the Ossenkopf et al. (1992) data, but they pro- 
duce an excess of flux at smaller wavelengths for a grain ra- 
dius of 0.1 /im. From the investigation of the effects resulting 
from a variation of Ogi in Sect. 4.2.3, we know that the flux 
at short wavelengths decreases with increasing grain radius, 
and that the decline of the visibility becomes steeper. There- 
fore, we have calculated a series of models with the David & 



Pegourie (1995 1 data where we varied the grain radius from 
flgr = 0.1 /im to 0.5 /im. Figure A. 3 shows the calculated SED 



and the 2.11 /im visibilities. The derived properties of the cor- 
responding models are given in Table |A.2| . Figure A.. 4 shows 
the extinction coefficient per unit volume of the grains Kcxt / V^r 
for the different grain radii. 

Compared to the models calculated with the sil-Ow con- 
stants, the sil-DP models show a flat flux distribution below 
A ^ 8 /tm with a steep drop at a certain wavelength, which 
depends on the grain radius. Furthermore, much higher values 
of Ogr are needed to match the observation. The short wave- 
length photometry of AFGL 2290 constrains the grain radius 
to flgr ^ 0.5 /tm, although the slope of the observed SED at 
wavelengths ^ 2.2 /im is less well reproduced compared to 
the sil-Ow model with agi =0.1 /im. The different behaviour 
below A = 8 /im is caused by the fact that the sil-DP con- 
stants yield smaller extinction efficiencies for grains of equal 
radius compared to the sil-Ow constants. Therefore, substan- 
tially larger grains are needed to produce comparable optical 
depths at shorter wavelengths. As a consequence, the contri- 
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Fig. A.l. SED (a) and visibilities (b) for models with the pa- 
rameters of model A, but with different dust optical proper- 
ties: 'warm' silicates from Ossenkopf et al. (1992) (solid line). 



'cold' silicates from Ossenkopf et al. (1992 ) (long dashed line) 



Draine & Lee (1984 ) (short dashed line), and David & Pegourie 
(1995) (dotted line). The corresponding optical depths and de- 
rived model properties are given in Table A.l The spectra 
have been scaled with different /b to match the observations 
at A > 8 ^m. 



bution of scattering to the extinction is still important up to 
wavelengths, which are approximately 3-4 times larger. 

Due to the stronger dependence of Kcxt / V^i on the grain ra- 
dius (Fig. A.4) and the resulting larger variation of the 2.11 fj,m. 



optical depth (see Table A.. 2), the changes of the 2.11 /itm visi- 
bility are much more pronounced. With increasing Cgr the cur- 
vature of the visibility increases, and V2.11 falls off to smaller 
values at q = 13.5arcsec~^ because the stellar contribution 
to the monochromatic flux is reduced. Up to Cgr « 0.2 fim 
the change of the curvature is mainly caused by the changing 
slope of the spatial intensity distribution, as discussed in Sect. 
4.2.3 . At larger grain radii the increase of the radius of the in- 
ner boundary and the corresponding reduction of the spatial 
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Fig. A.2. Extinction coefficient per unit volume of the grains 
Kcxt/V^r for flgr = 0.1 /im obtained from different opti- 
cal constants: 'warm' silicates from Ossenkopf et al. ( 1992| ) 



(solid line), 'cold' silicates from Ossenkopf et al. (1992) (long 
dashed line), Draine & Lee (1984i (short dashed line), David 



& Pegourie (1995 1 (dotted Une) 



Table A.2. Derived properties for mo dels with the optical con- 
stants from David & Pegourie ( 1995 ) and different grain radii. 
The other parameters are equal to model A, except for ro.55 
which has been adjusted to fit the observed SED. 
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frequencies, where V2.11 approaches a constant value (or zero) 
becomes more important. 

The best match of the observed visibility V^bs is obtained 
for the model with Ogr — 0.2 /im (dotted line in Fig. A. 3), 
which is somewhat larger than a. 



= 0.16 /im for the best 
matching model with the sil-Ow optical constants. The result- 
ing dust mass loss rate is about 30 % smaller, and the radius of 
the inner boundary is about 12 % smaller. Again the curvature 
of the model visibility is slightly too strong. But in contrast to 
the model with the sil-Ow optical constants, there is now an 
excess of flux at smaller wavelengths A ^ 5 /im. 



As discussed at the end of Sect. 1-.2.3, we would expect a 
deficiency of flux if the CDS of AFGL 2290 has a disk-like 
structure. This suggests that the optical constants from David 



& Pegourie (1995) underestimate the extinction of the grain 
material at short wavelengths. However, we also do not k now, 
whether the optical constants from Ossenkopf et al. ( 1992 ) rep- 
resent the intrinsic optical properties of the grains, because the 
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Fig. A.4. Extinction coefficient per unit volume of the grains 
Kcxt/V^i calculated with the optical constants from David & 
Pegourie ( |1995 ) for different grain radii ranging from Ogr = 
0.1 /im to flgr — 0.5 /im. 
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Fig. A.3. SED (a) and visibilities (b) for models with the opti- 
cal constants from David & Pegourie (1995 ) and different grain 
radii ranging from Cgi = 0.1 /im to flgr = 0.5 /^m. The other 
parameters are equal to model A. The corresponding optical 



depths and derived model properties are given in Table A. 2. 
The spectra have been scaled with different /b to match the 
observations at A > 8 ^m. 



shape of the SED at short wavelengths is affected to an un- 
known degree by the geometry of the non-spherical dust dis- 
tribution around AFGL 2290. For the same reason we refrain 
from enforcing a match of the SED at short wavelengths by a 
suitable modification of the dust optical constants in this re- 
gion. 



